Pairing Gap and Polarisation Effects 
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The phenomenological adjustment of the nuclear pairing strength is usually performed with respect 
to the odd-even staggering of the binding energies. We find that the results strongly depend on the 
way in which the ground states of the odd nuclei are computed. A thorough calculation including 
all time-even and time-odd polarisation effects induced by the odd nucleon produces about 30% 
reduced odd-even staggering as compared to the standard spherical calculations in the relativistic 
mean-field model. The pairing strength must be enhanced by about 20% to compensate for that 
effect. The enhanced strength has dramatic consequences for the predicted deformation properties 
of the underlying mean-field models, possibly implying that new adjustments of their parameters 
become necessary as well. 

PACS numbers: 21.30.Fe, 21.60.-n, 24.10.Jv 



Pairing is an essential ingredient of nuclear mean-field 
models since the early days of the single-particle shell 
model [0,0| . It can be motivated from a theoretical point 
of view by a near-singularity of the effective T matrix 
near the Fermi surface 0] and it manages to subsum- 
marize a large amount of two-body correlations at the 
price of a moderately extended mean-field theory. From 
a phenomenological point of view, pairing is needed to 
explain, e.g., the existence of spherical non-magic nuclei 
and to account for the observed strong odd-even stag- 
gering of the binding energies. The purely theoretical 
access to pairing from nuclear many-body theory is still 
plagued by unresolved quantitative problems of nuclear 
"ab initio" calculations H . One thus has to recur to phe- 
nomenological information which is exploited nowadays 
with a high degree of systematics (for a recent and com- 
prehensive compilation see H). The general strategy of 
these phenomenological evaluations is to relate the odd- 
even staggering of the binding energies to the pairing gap 
by taking appropriate differences, e.g., for neutrons the 
fourth order difference 

A(f) = -|(-1)^ [E(N - 2) - 4.E{N - 1) + &E{N) 

-AE{N+l) + E{N + 2)\ (1) 

at fixed proton number. This access, however, has the 
basic problem that pairing is not the only source of odd- 
even staggering. Odd-even fluctuations in energy can also 
be produced by polarisation of the even core through the 
odd nucleon, which leads to odd-even jumps in deforma- 
tion, spin- alignment or dynamical currents. For example, 
there is pronounced odd-even staggering in the energies 
of metal clusters which, however, can be explained ex- 
clusively by effects of spin-alignment and Jahn- Teller de- 
formation B. Spin alignment can be excluded for most 
nuclei. The Jahn- Teller effect alone can produce size- 
able odd-even effects in deformed nuclei as was worked 



out in a recent paper H], but its influence is diminished 
for semi-magic nuclei if pairing is switched on because 
this restores spherical symmetry over the whole isotopic 
or isotonic chain. Nonetheless, there remains a polari- 
sation of the core through the field of the odd nucleon, 
which may deliver substantial contributions to the odd- 
even staggering and thus should be taken into account 
when adjusting pairing strengths to the phenomenologi- 
cal pairing gap (|l|) . The question is how large these polar- 
isation effects are in practice. There is, first, the static 
deformation-polarisation which is induced by the finite 
multipole moment of the odd nucleon's density. In ad- 
dition the odd nucleon breaks the intrinsic time-reversal 
invariance because it carries a nonzero spin and a current 
contribution. These time-odd components can induce a 
sizeable time-odd response in terms of spin and current 
polarisation in the even-even core |q] . In the following we 
call that a dynamical polarisation to distinguish it from 
the mere deformation effects. It is the aim of this letter 
to investigate the effect of such polarisation effects on the 
odd-even staggering and thus on the adjustment of the 
pairing gap (|^). 

The investigation requires a mean-field model which 
provides reliable response properties in all channels, in- 
cluding unnatural parity states and spin polarisation. 
The relativistic mean-field (RMF) model includes a "nat- 
ural" description of spin properties |p|-pO[ and is a reliable 
starting point for our present investigation of dynami- 
cal polarisation. A similar investigation on the grounds 
of the non-relativistic Skyrme-Hartree-Fock model was 
performed recently fO]. We employ the RMF with the 
parametrization PL-40 p2[ similar as in a previous sys- 
tematic investigation of odd nuclei near magic shells [^ . 
The details of the model, the strategy to define the block- 
ing in odd nuclei, and the numerical handling can be 
found there. We perform axially symmetric deformed 
calculations including time-odd currents to take into ac- 
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-346.0 
-346.0 
-346.0 
-412.4 



1.2587 
1.2493 
1.1600 
0.8104 
1.2588 



TABLE I. Neutron gaps Ai"*' in ^^"^Sn for the various levels 
of approximation. In case of dynamical polarisation the gaps 
are shown for the original as well as for the refitted pairing 
strength V„ (last line). 



count the broken time-reversal symmetry and its related 
polarisation effects. Pairing is described at the BCS level 
using matrix elements computed from a zero-range pair- 
ing force Vq 6{ri — r2) with q G {p, n} acting in a pairing 
space which is limited by a soft cut-off of Woods-Saxon 
shape in energy [ ^3^1^ ] with the cut-off adjusted dy- 
namically to include Nq + 1.65 Nq' nucleon states [Oj. 
The pairing strengths are fitted to the experimental gaps 
within several semi-magic isotope and isotone chains in 
connection with the actual mean-field force. The odd nu- 
clei needed for this approach are calculated in a spherical 
approximation. For PL-40, which is used in this study, 
we find Vp = -348 MeVfm^ and K = -346 MeVfm^. 
These values correspond to a standard determination of 
pairing strengths without polarisation effects. We take 
them as a base point for comparison with the more elab- 
orate adjustments to be discussed in the following. 

As a test case we consider the neutron gaps in the chain 
of Sn isotopes. The magic proton number in Sn allows to 
concentrate on the neutron gaps and it renders all even 
isotopes spherical, which confines the deformation effects 
exclusively to the odd isotopes. 

Figure shows the results for the neutron gap (nj) at 
various levels of mean- field calculations: spherically, de- 
formed but time-even, and deformed with dynamical po- 
larisation (i.e., including time-odd currents). As men- 
tioned above, the spherical results had been fitted to 
agree with the experimental gaps in the average, which 
has been achieved more or less nicely. Allowing for de- 
formation (yet without dynamical polarisation) indeed 
changes the gaps at several isotopes. For the heavier 
nuclei, however, to which pairing is usually fitted, the 
resulting reduction of the gap is sufficiently small to ne- 
glect them at the level of quality with which pairing can 
be adjusted anyway. This can be deduced from the small 
changes from "spherical" to "deformed" seen in fig. n. 
The effect is much larger when allowing for dynamical 
polarisation. The gaps shrink in average by about 30%. 
The example here shows that dynamical polarisation can 
have a large effect on the odd-even staggering which, in 
turn, puts a warning signs on a phenomenological ad- 
justment of pairing properties on the the basis of energy 
differences. One should, in fact, adjust the gaps (|^) anew 
while using fully polarised calculations. 
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FIG. 1. Fourth-difference neutron gaps (hi) for Sn isotopes 
{Z = 50) computed with the RMF parametrization PL-40 in 
different approaches. Spheres stand for purely spherical cal- 
culations without time-odd currents, triangles for deformed 
calculations without time-odd currents, and squares for un- 
restricted deformed calculations including time-odd currents. 
Experimental gaps are given as full rhombi for comparison 
(the masses were taken from [[l6|). 



In order to check that this reduction of the gap is not 
a particular feature of PL-40, we have performed simi- 
lar calculations for the parametrization NL3 fl^ which 
was adjusted with different bias, we find for the case of 
^^^Sn the sequence 1.32 MeV for the gap from spherical 
calculations, over 1.30 MeV for the gap from deformed 
and time-even calculations down to 0.93 MeV for the case 
with full dynamical polarisation. That is very much the 
same reduction as explored for PL-40, see Table |. we 
thus are confident that the effect is of generic nature, at 
least for the RMF. 

We now want to exemplify the consequences of such a 
modified pairing strength for the present test case. To 
that end, we have performed a readjustment of the pair- 
ing strength taking properly into account all polarisation 
effects. In order to cope with the rather large expense 
of the full-fledged calculations of odd nuclei, we have de- 
cided to concentrate on the neutron pairing strength and 
to refit that with respect to one nucleus, namely ^^^Sn. 
We find that the readjustment (including time-odd cur- 
rents) increases the required pairing-strength parameter 
by as much as 19%. Proton pairing will be needed in 
two applications later on. We find that the same en- 
hancement factor of 19% reproduces gaps in heavy nu- 
clei very well. For that we assume the same upscaling 
factor as was found for neutron pairing. This minimal 
fitting strategy is sufficient for the present purposes of 
an exploratory study. The results for the case of neutron 
pairing in ^^^Sn are summarised in table |. One sees 
that the dynamical polarisation reduces the gap by 35% 
in this case and a counteracting readjustment increases 
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FIG. 2. Ground state observables in Sn isotopes drawn 
as functions of the neutron number N. Upper panel 
two-neutron separation energies 82,1 = E{N — 2) — E{N) 
second from above: rms radii of the charge distribution, R'^^^ 
third panel: diffraction radii R^^; lowest panel: surface 
thickness a. The error bars represent the average precision of 
the mean field model to reproduce these quantities (the ac- 
tual experimental errors would be much smaller). Two cases 
are compared: a calculation with the spherically fitted pair- 
ing strength (full lines) and with the refitted strength (dashed 
lines). 



the pairing-strength parameter by 19%. It remains to 
check the consequences for other observables. 

A summary of ground-state properties is shown in 
fig. g. The binding energies are displayed in terms 
of the two- neutron separation energies S2niN,Z) = 
E{N — 2,Z) — E{N, Z) because this amplifies possible 
effects. The overall size and trend is not much affected 
by the change of the pairing strength. There is, however, 
a systematic modification near the shell closures. Con- 
sider, e.g., the magic N = 82. The larger, refitted pairing 
strength produces a smaller S2n before this shell and a 
larger S2n right after it. The jump of the two-neutron 
separation energies at a magic shell serves as an empir- 
ical measure of the "magicity" , called the two-neutron 
shell gap. And we see that this shell gap is reduced 
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FIG. 3. Quadrupole deformation properties for ^^"Sm 
drawn versus the dimensionless quadrupole deformation 
/32 = [47r/(3Ai?2)] (r^Yao) with R = 1.2 fm A^''^ . The full 
line shows the result obtained with standard pairing (fitted to 
spherically calculated nuclei only) and the dashed line those 
for the enhanced pairing strength (fitted including polarisa- 
tion effects). Upper: Deformation energy, i.e. total energy 
rescaled to zero at the spherical minimum. Middle: Neutron 
pairing energy. Lower: Proton pairing energy. 



by as much as 25% when employing the refitted pair- 
ing strength. The binding energies as such are slightly 
enhanced by at most 0.25% in the mid-shell region due 
to the stronger pairing. This change can be, and needs 
to be, corrected by a slight readjustment of the RMF 
parametrization. 

Fig. pi also shows the bulk properties of the nuclear 
charge distribution, i.e., root-mean-square radii [nSJ, 
diffraction radii, and surface thicknesses |13]. There are 
small effects visible but they remain far below the pre- 
cision of the force to describe these observables (indi- 
cated by the error bars). Density and formfactor are thus 
rather robust against this change in pairing strength. 

The situation may be different for deformation proper- 
ties because these are known to result from an interplay 
between shell structure and pairing. As an example, in 
fig. y we show the quadrupole deformation energy for 
^^°Sm, the softest member of the Sn isotope chain. The 
results confirm the rule-of-thumb that stronger pairing 
acts more to restore the spherical shape. The refitted, 
enhanced pairing strength clearly produces a stiffer de- 
formation energy curve and suppresses more efficiently 
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FIG. 4. Deformation energy curve (defined as in as fig. y) 
for ^''"Pu. The points with error bars show the experimental 
minima and barriers [GOli 



the deformed side minima which are visible in case of 
the standard pairing strengths. It is worthwhile to esti- 
mate the consequences for the low-lying 2+ state. The 
curvature at /?2 = is enhanced by about 130%. The 
cranking mass is proportional to the inverse quasiparti- 
cle energies which are enhanced by about 35% for the 
refitted pairing strengths. The mass is thus reduced by 
25% and the energy of the low-lying 2+ state is then al- 
together increased by about 60% for the refitted pairing 
strengths. This little estimate demonstrates that an ap- 
propriate determination of the pairing strength is crucial 
for these low-lying collective modes. 

The two lower plots of fig. g show the proton and neu- 
tron pairing energies versus deformation. The larger re- 
fitted strength yields, of course, much larger pairing en- 
ergies (absolute values), but it does also produce a larger 
diffuseness of occupation numbers near the Fermi surface 
which, in turn, degrades the mean-field binding leaving 
altogether the net effect of a stiffer deformation energy 
curve as seen in the upper panel of the figure. The ex- 
ample does also show the breakdown of pairing towards 
the proton shell closure at spherical shape, and it shows 
that the breakdown comes, of course, a bit later for the 
stronger pairing. This feature, however, should not be 
overstressed because the breakdown is anyway an arte- 
fact of the mere BCS treatment. A more correct ap- 
proach would be to use particle number projection. But 
the main effect on the deformation energy curve will not 
be changed by that. 

Other important observables in nuclear physics are fis- 
sion barriers and it is worth looking at the effect when 
the pairing strengths are changed. Fig. shows the en- 
ergy curve for asymmetric fission of ^^'^Pu. For this heavy 
nucleus, we see even more dramatic changes. The (de- 
formed) minima are softened such that the curvature is 
reduced by a factor of two. At first glance, one may 
be surprised to see here a softening whereas the higher 



pairing strength caused a stronger curvature in case of 
^^"Sn. It is to be noted that the minima in these two cases 
have different origins. For ^^°Sn, the high level density 
at spherical shape inhibits, in principle, a spherical mini- 
mum and sphericity is only restored by the action of pair- 
ing. Consequently, stronger pairing makes the spherical 
shape even more pronounced. For ^'^''Pu, on the other 
hand, the deformed minimum is caused by a low level 
density at this place and pairing can only counteract this 
preference, thus delivering a softening of the minimum. 
There may be, nonetheless, only a small effect on the 
energy of the vibrational states in the deformed minima 
because the quasiparticle energies (and with them the 
Inglis mass) are enhanced also by about a factor of two. 
On the other hand there is a strong effect on fission prop- 
erties because the barriers are lowered by 1-2 MeV. The 
change of the barriers moves them away from the experi- 
mental points for this particular mean-field parametriza- 
tion (note that the inner barrier would be further lowered 
when allowing for triaxial shapes). This means that the 
selection of the appropriate parametrization for fission 
p2| needs to be reconsidered in connection with the new 
pairing strength. 

Altogether, we find that dynamical polarisation effects 
can have a strong influence on the pairing gap as de- 
duced from even-odd staggering of binding energies. In 
the present test case, one needs to enhance the under- 
lying pairing strengths by about 20% to compensate for 
that effect. This change in pairing strengths has only a 
small effect on the ground state properties of even-even 
nuclei like binding energies and radii. But more elabo- 
rate quantities like the jump of the two-particle separa- 
tion energies at magic shells can react sensitively. Even 
more dramatic changes are seen for deformation proper- 
ties. Vibrations around spherical shapes become more 
rigid and fission barriers are significantly lowered. These 
findings are disquieting and call for further critical in- 
spection with different test cases and other mean-field 
theories. There is, for example, one possible fiaw in the 
RMF. It is an effective Hartree theory. Thus the po- 
larisation energy in the odd nucleus contains a contri- 
bution from the self-interaction of the odd nucleon and 
a self-interaction corrected theory may produce different 
results. The situation is different in the non-relativistic 
Skyrme-Hartree-Fock model which can be considered as 
a true Hartree-Fock variational theory and which thus 
is free from the self-interaction effect. Research in this 
direction is underway. 
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